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ABSTRACT 

The sub-arcsec bright points (BP) associated with the smah scale magnetic fields in 
the lower solar atmosphere are advected by the evolution of the photospheric granules. 
We measure various quantities related to the horizontal motions of the BPs observed in 
two wavelengths, including the velocity auto-correlation function. A 1 hr time sequence 
of wideband Ha observations conducted at the Swedish 1-m Solar Telescope {SST), and 
a 4 hr Hinode G-band time sequence observed with the Solar Optical telescope are used 
in this work. We follow 97 SST and 212 Hinode BPs with 3800 and 1950 individual 
velocity measurements respectively. For its high cadence of 5 s as compared to 30 s 
for Hinode data, we emphasize more on the results from SST data. The BP positional 
uncertainty achieved by SST is as low as 3 km. The position errors contribute 0.75 
km^ s~^ to the variance of the observed velocities. The raw and corrected velocity mea- 
surements in both directions, i.e., {vx,Vy), have Gaussian distributions with standard 
deviations of (1.32,1.22) and (1.00,0.86) km s~^ respectively. The BP motions have 
correlation times of about 22 — 30 s. We construct the power spectrum of the horizontal 
motions as a function of frequency, a quantity that is useful and relevant to the studies 
of generation of Alfven waves. Photospheric turbulent diffusion at time scales less than 
200 s is found to satisfy a power law with an index of 1.59. 



Subject headings: Sun: photosphere — Sun: surface magnetism 
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Introduction 



The discrete and small scale component of the solar magnetic f ield is revealed in the high spatial 



resol ution observations of the Sun. Ground based observations (IMulleij Il983l . Il985l : iBerger et al 



19951 ) show clusters or a network of many bright points (hereafter BPs) in the inter-granular lanes, 
with each individual BP having a typical size of 100 — 150 km. These BPs are known to be kilogauss 
flux tubes in the small scale magneti c field (SMF) , and are extensively used as proxies for such 



flux tubes dC hapma n &: SheelevI Il968l : IStenflo 



1973 



Stenfio fc Harvev 



1985 



Title et al. 



19871 . see 



de Wijn et al., (.2009. ) for a review on the SMF). High cadence observations and studies show that 
magnetic BPs are highly dynamic and intermittent in nature, randomly moving in the dark inter- 
granular lane^. These motions are mainly due to the buffeting of granules. The SMF is passively 
advected to the boundaries of supergranules creating the magnetic network in the photosphere. 

Earlier works by several authors have reported mean rms velocities of magnetic elements in 
the or der of a few km s~^ . With the ground based observations of the granules at 5750 A (white 
light), iMuller et al.l (119941) have identified many network BPs with turbulent proper motion and a 



mean speed of 1.4 km s . iBerger &: Titld (119961 ) have used G-band observations of the photosphere 



and found that the (j- band BPs move in the intergranular lanes at speeds from 0.5 to 5 km 



Berger et al.l (|l998l ) observed the flowfield properties of the photosphere by comparing the 
magnetic network and nonmagnetic quiet Sun. They show that the convective flow structures 
are smaller and much more chaotic in magnetic region, with a mean sp eed of 1.47 km s~^ fo r 



the tr acked magnetic BPs. With the G-band and continuum filtergrams, Ivan Ballegooijen et al 



(|l998l ) used an object tracking technique and determined the autocorrelation function describing the 
temporal variation of the bright point veloci ty, with a correlation time of about 100 s. Correcting for 
measurements errors. iNisenson et al.l ( 20031 ) measured a 0.89 km s^^ rms velocity of BPs. Advances 
in the ground based observations like rapid high cadence sequences with improved adaptive optics 
to minimize seeing effects, and also the space based observation s at high resolut ions, continued to 
attract many authors to pursue BP motion studies. For example. lUtz et al.l (120101 ) used space based 
Hinode G-band images to measure BP velocities and their lifetime. The BP motions can be used 
to measure dynamic properties of magnetic flux tubes and their interaction with granular plasma. 
Photospheric turbulent diff usion is one such dynamic al aspect that can be derived consequently 
from the BP random walk. iManso Sainz et al.l ( 2011) measured a di f fusion constant of 195 km^ 
s~^ from the BP random walk and their dispersion. lAbramenko et al.l ( 201ll ) studied photospheric 
diffusion at a cadence of 10 s with high resolution TiO observations of a quiet Sun area. They 
found a super-diffusion regime, satisfying a power law of diffusion with an index 7 = 1.53, which is 
pronounced in the time intervals 10 — 300 s. 



The implications of these magnetic random walk motions have been found very fruitful recently. 



^See an accompanying movie of SST data used in this study. Also see 

TiO movies at |http : / /www . bbso . n j it . edu/ iist_gallery . html [ and G-band movies at 
jhttpT7/solar-b . nao . ac . jp/qLmovies/index_e . shtml . | 



- 3 - 



Such motions are capable of launching magneto- hydrodynamic (MHD) waves (|Spruitlll98ll ). which 
are potential candidates for explaining the high temperatures observed in the solar chromosphere 



and c orona. For example, a three-dimensional MHD model developed by Ivan Ballegooiien et al 



(|201ll ) suggests that random motions insid e BPs can create Alfven wave turbuleii ce, which dissi- 
pa tes the waves in a coronal loop (also see lAsgari-Targhi &: vaii BallegooiienI |2012| ) . Observations 



bv lDe Pontieu et alJ (j2007b|) ; |jess_etalJ ([20091); 



Mcintosh et al 



(j201ll ) provide strong evidence that 



the Alfvenic waves (which are probably generated by the BP motions), have sufficient energy to 
heat the quiet solar corona. To test theories of chromospheric and coronal heating, more precise 
measurements of the velocities and power spectra of BP motions are needed. 



Nisenson et al.l (|2003l ) worked on the precise measurements of BP positions, taking into account 
the measurement errors. The auto-correlations derived by them for the x— and y— components 
of the BP velocity using high spatial resolution and moderate cadence of 30 s observations gave 
a correlation time of about 60 s, which is twice the cadence of the observations. This suggests 
an over-estimation of correlation time and an under-estimation of the rms velocity power, with a 
significant hidden power in the time scales less than 30 s, and thus warranting for observations at 
even higher cadence. This is important because, the measured power profile, which is the Fourier 
transform of the auto-correlation function, gives us an estimation of the velocity amplitudes and 
energy flux carried by the waves that are generated by the BP motions in various and especially at 
higher frequencies. 

In this study, we use 5 s cadence wideband Ha observations from the Swedish 1-m Solar 
Telescope to track the BPs and measure their rms velocities. For comparison, we also use a 30 
s cadence G-band observational sequence from Solar Optical Telescope onboard Hinode. These 
independent and complementary results take us closer to what could be the true rms velocity and 
power profile of the lateral motions of the BPs. The details of the datasets used, analysis procedure, 
results and their implications are discussed in the following sections. 



2. Datasets 



In this study, we have analyzed time sequence of intensity filtergrams with 5 and 30 s cadence. 
A brief description of the observations is given below. 

5 sec data: These observations were obtained on 18 June 2006, with the Swedish 1-m ^o- 
lar Telescope (^ST; Scharmer et al. 2003a) on La Palma, using the adaptive optics system (AO; 
Scharmer et al. II2OO3I) in combination with the Multi-Object Multi-Frame Blind Deconvolution (MOMFBD; 



van Noort et al 



20051 ) image restoration method under excellent seeing conditions. The target area 
is a quiet Sun region away from disk center at {x,y) = (—307", —54") and fi = 0.94 (see Figured]). 
The time sequence is of one hour duration starting at 13:10 UT. Here we analyze images froin 
the wideband channel of the Solar Optical Universal Polarimeter (SOUP; Title Sz Rosenberg 1981 ) 
which received 10% of the light before the SOUP tunable filter but after the SOUP prefilter (see 
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De Pontieu et al.ll2007al . for the optical setup of the instrument). The prefilter was an FWHM=8 
A wide interference filter centered on the Ha line. The SOUP filter was tuned to the blue wing 
of Ha at -450 mA but that data is not considered here. On the wideband channel there were 
2 cameras (running at 37 frames per sec) positioned as phase-diversity pair - one in focus and 
one camera 13.5 mm out of focus, the data from the two cameras has been processed with the 
MOMFBD restoration method in sets of 5 seconds, creating a 5 s cadence time sequence with a 
total of 720 images. After MOMFBD processing, the restored images were de-rotated to account 
for the field rotation due to the altazimuth mount of the telescope. Furthermore, the images were 
aligned using cross-correlation on a large area of the field of view (FOV) . The images were then 
clipped to 833 x 821 pixels (with 0.065" per pixel), to keep the common FOV (the CCDs have 
1024 X 1024 pixels, some pixels are lost after alignment between focus and defocus cameras). 

For a reference direction, the solar north in the SST time sequence is found by aligning an 
earlier SST observation of that day of a magnetogram of an active region (AR) to a full disk 
SOHO /MDI magnetogram (the AR was just outside the MDI High Res region). From that com- 
parison, we fix the direction of solar north and disk center (black and white arrows respectively in 
Figured]). Though we do not rotate the images to match the solar north during our analysis, the 
angles are taken into account at a later stage to correct for the projection effects in the velocity 
measurements. 

30 s ec data: We use (?-band filtergranis obse rved with Solar Optical Telescope (SOT) on board 
Hinode (jKosugi et al.l 120071 : iTsuneta et al.l l2008l ) , on 14 April 2007. The observations were made 
for a duration of 4 hr, with a 30 s cadence in a FOV of 55" x 55" ( 0.05" per pixel; 1024 pixels in 
both X and y directions), near disk center. The images were processed using standard procedures 
available in the solarsoft library. 



3. Procedure 



In this section, we briefiy describe the method of determining the BP positions, and the velocity 
measurements through the correlation tracking. 



3.1. BP Positions 

We manually select the BPs to estimate their position to a sub-pixel accuracy. We consider 
the coordinates of maximum intensity of a given BP to be the position of that BP and the method 
for measuring these positions involves two steps. In the first step, we visually identify a BP and it 
is selected for analysis for a period during which it is clearly distinguishable from the surrounding 
granules. On an average, we follow a BP for about 3 — 5 min. The BPs with elongated shapes 
are not considered for the analysis. Also, we stop following a BP if it is substantially distorted or 
elongated from its initial shape. Though time consuming, manual selection gives a handle on the 
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validity of the positional accuracy of a BP from frame to frame. At each time step, using a cursor, 
an approximate location ix'app,yapp) of a particular BP is fed to an automated procedure to get its 
accurate position, which is step two in our method. 

Step two is completely an automated procedure. Here, we use a surface interpolation technique 
to get a precise position of that BP (to a sub-pixel accuracy). Approximate position from the 
previous step is used to construct a grid of 5 x 5 pixels covering the full BP (with (xjjpp, y^^p) as 
the center of that grid). Now, our procedure fits a 2-D, 4*^^ degree surface polynomial to that grid 
(using SFIT, an IDL procedure); interpolates the fit to one- hundredth of a pixel; returns the fine 
location of its peak {6x' , 5y') within that grid and finally stores the accurate position {x'^p^y'^p) of 
that BP (which is the sum of its approximate and fine positions (a^^pp + 5x' , y'app + ^u')), for further 
analysis. Therefore, the position of a BP with index j in a frame i is given by 

(x'bp, y'Bp)i = (x'app + Sx',y'^pp + 5y')l (1) 
and all the coordinates till this point are relative to the lower left corner of the image. 



3.2. Reference Frame 

Though the positional measurements of BPs as described in Section 13.11 are accurate, they 
cannot be directly used to measure the velocities as there are artificial velocity sources viz. the 
instrumental drifts, seeing variations, jittery motions and also the solar rotation, which vectorially 
add to BP velocities and thus are required to be removed from the analysis. While Hinode (space 
based) data is not subjected to seeing variations, SST (ground based) data has been corrected for 
seeing as described in Section [2j Further, we need to correct for instrumental drifts, jitters and 
solar rotation. Calculating the offsets between the successive images is necessary to remove these 
artificial velocities. In this section, we describe the method of our cross correlation analysis used 
to co-align the images. 

Cross correlation (C) of two images f{x,y) and g{x,y) is defined as 

^ ^ 1 ifix,y) -J)ig{x,y) - g) 

x,y J y 

where, / \g] and cj/ [ug] are mean value and standard deviation of f(x,y) [g(x,y)] respectively, 
k is the number of pixels in each image, for normalization. With the above definition of cross 
correlation, to get the offsets between two images, we need to shift one image with respect to the 
other (in both x and y directions ) and find at what offsets (independent in x and y) the correlation 
function attains the maximum value. In general, for shifts of — / to +1, the cross correlation is a 
2-D function with 21 + 1 rows and columns. Let Ix and ly be the coarse offsets between the two 
images in x and y directions, respectively, such that the cross correlation reaches its maximum 
value: inax{C) = C{lx,ly), where —I < lx,ly < I- To get the sub-pixel offsets, the fine offsets 
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{5lx, Sly) are calculated. The method is similar to finding the fine position of BP by using a 5 x 5 
pixel grid but now about {lx,ly) of C. 

Instead of cross correlating every successive image with its previous one, we keep a reference 
image for about 200 s, i.e., a frame i taken at time t (it) is used as a reference for the subsequent 
frames till i+200 s {it+200) for cross correlation. Therefore the 5 sec (SST) and the 30 sec ( Hinode) 
data have about 40 and 7 images respectively in each set. By keeping the last image of a set equal 
to the first image in its next set, we can co-align different sets. In this way, the accumulation of 
errors in the offsets can be minimized. 

With the above background on co-aligning images to find various drifts, we present the results 
of drifts found in Hinode data. As an illustration, we divide the full (i.e., 55" x 55" x 4 hr) Hinode 
time sequence into 4 quadrants with 27.5" x 27.5" x 4 hr each. Further, we do correlation tracking 
(as described above by keeping 7 frames per set) on each quadrant separately and plot the results 
in Figure [2j The four dashed lines in the left and the right panels are the offsets in x— and y— 
directions respectively, the thick dashed line in each panel is the average of the offsets (i.e., average 
of four dashed lines), and the solid red curve is the offset obtained by considering the full FOV. 
Clearly, in each quadrant, the offsets have a trend similar to that of the full FOV (solid red curve) 
and an additional component of their own. This additional component is probably the real velocity 
on the Sun due to flows with varying length scales (for example, super-granular, meso-granular and 
granular) and with flow directions changing over areas of a few tens of arcsec^ on the Sun. 

In this paper we are mainly interested in the dynamics of the BPs relative to their local 
surroundings, as granulation flows will have a dominant effect on the BP velocities and their 
variations on short time scales. Hence, we consider a 5" x 5" area about the BP as a reference 
frame for that BP (i.e., keeping the BP in the center of the local area). The cross correlation is 
performed on this 5" x 5" area instead of the full FOV to get the offsets, which are subtracted from 
the {x'^p,y'Bp)i- The BP positions corrected for offsets are now given by 



where jiocai represents the local area of BP-' . 

In the case of the SST data, the observations are off disk center at (—307", —54"), which 
corresponds to a heliocentric angle of arccos(0.94). This will introduce a projection effect on the 
measured horizontal velocities in both x'— and y'— directions and needs to be corrected. To do 
this, the coordinate system {x',y') defined by the original SST observations, is rotated by 45° in 
the anti-clockwise direction. Now the image plane is oriented in E-W (parallel to equator, new 
X—) and N-S (new y—) directions. Further, the E-W coordinate is multiplied by a factor of 0.94~^. 
Hence the new coordinate system (x, y) is given by 



''BPC^y'BPc)l — i^BP^ UBpYi ~ i^x + Slx,ly + SlyY- 



local 



(3) 



X = (x' cos 45° -|- y' sin 45°) X 
y = (—x' sin 45° -|- y' cos 45°). 



1 



0.94' 



(4) 
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SST BP positions {x'^pc^ y'BPc)h measured from Equation ([3]), are remapped to {xbpc-, xbpc)1 
using the above coordinate transformation^. 



i 1 



4. 



Results 



In this section, we present various results in detail giving more emphasis on the SST results. 
We have selected 97 SST BPs with ~3800 individual velocity measurementjf]. Figure [3] shows 
the paths of four individual SST BPs. Some of the BPs move in a relatively smoother path 
while some exhibit very random motions to the shortest time steps available. BPs drift about 
a few hundred km in a few minutes. The instantaneous velocity {vxiVy)lj^i of a BP is given by 
{xbpc, yBPc)i+i — {xbpCi yBPcYi, multiplied by a factor to convert the units of measured velocity 
to km s~^ (9.4 in case of SST which is, image scale of SST in km divided by the time cadence in 
sec). Figure m shows the plot of such velocities as a function of time for BP^^S (path of BP#3 is 
shown in the lower left panel of Figure [3|) . Usually, the changes in the velocity are gradual in time 
but, sometimes we do see sudden and large changes in the magnitude and direction of the velocity 
(for example at 1 min in Vx and at 2 min in Vy in Figured]). Note that a large change of velocity of 
one sign is followed immediately by a change of the opposite sign, so the net change in position is 
not very large. This suggest that these changes are due to errors in the positional measurements. 
A position error at one time will affect the velocities in the intervals immediately before and after 
that time. In the following we will assume that such changes in velocity are due to measurement 
errors. However, we cannot rule out that some of these changes are due to real motions on the Sun 
on a time scales less than 5 seconds. 

The means and standard deviations of Vx and Vy are listed in Table [T] (first line). Histograms 



of the distribution of velocities Vx-, Vy and v = yu^ + are shown in Figure [5] (panels (a), (b) 
and (c) respectively). Solid lines in panels (a) and (b) are Gaussian fits to the histograms with raw 
standard deviations (cr^^r) of 1.32 and 1.22 km s^^. A scatter plot of Vx against Vy is shown in panel 
(d), which is symmetric in the f— space. However, a small non-zero and positive mean velocity of 
about 0.2 km s~^ is noticed, suggesting that there is a net BP velocity with respect to the 5 arcsec 
boxes that we used as reference frames. Values of the mean and rms velocities as determined from 
the fits are also listed in Table [l] (second line) . These distributions are a mix of both true velocities 
and measurement errors. 

We can gain more insight into the the dynamical aspects of the BP motions by studying their 



^Note that the transformations in Equation (U) are only to modify the SST BP positions and in the rest of the 
paper, we use {x,y) for the remapped {x',y') of SST and {x,y) of Hinode. 

^Similarly, we have identified 212 Hinode BPs with 1950 individual velocity measurements. 
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observed velocity correlation function c{t), defined as 

Cxy.n = {vl^iVl^i+n)^ (6) 

where Cxx,m ^yy,n s-^'s the auto-correlations, and Cxy^n is the cross-correlation of Vx and Vy, n is the 
index of the delay time. (...) denotes the average over all values of the time index i and BP index 
i but for a fixed value of n. These results are shown in Figure [6l Top left and right panels are the 
plots of Cxx and Cyy respectively. Black curves are for the SST whereas the red curves show the 
Hinode results for comparison. Both the SST and Hinode results are consistent for delay times < 1 
min. However, the Hinode auto-correlations quickly fall to lower values. This is mainly a statistical 
error, since we do not have a large number of measurements in the case of Hinode. Focusing on 
periods < 1 min, it is clear from the auto-correlation plots that the core of the Hinode data within 
±30 s delay time, which is sampled with three data points is now well resolved with the aid of the 
SST data. Also, at shorter times, c takes a cusp-like profile. Extrapolating this to delay times of 
the order of 1 s, we expect to see a steep increase in the rms velocitiej^ of the BP motions. The 
bottom left panel shows the cross-correlation as function of delay time. The SST data show a small 
but a consistent and overall negative Cxy while the Hinode data show a small positive correlation. 
We suggest that the real cross-correlation Cxy ~ 0, and that the measured values are due to a small 
number of measurements with high velocities (largely exceeding the rms values). The lower right 
panel of Figure [6] shows the number of measurements A^^ used in the correlation analysis for both 
the SST and Hinode data. To obtain good statistics we collected enough BP measurements to 
ensure that Nn ^ 500 for all bins. 

In the rest of the section, we describe the method of estimating th e errors in the v elocit y 



measurements due to positional uncertainties by analyzing c{t). Following iNisenson et al.l (|2003l ). 
we assume that the errors in the positions are uncorrelated from frame to frame and randomly 
distributed with a standard deviation of ap. Since the velocities are computed by taking simple 
differences between position measurements (see above), the measurement errors increase the ob- 
served velocity correlation at n = by A (error), and reduce the correlations at n = ±1 by —A/2 
where A = 2{ap/6t)'^ and 6t is the cadence (see Equation 3 in their paper). We define 



A„ 



when n = 

iA when n = ±1 (7) 



2 



otherwise, 



which is valid only with our two-point formula for the velocity. Once A is determined, the rms 
values (crt,,c) of the true solar velocities can be measured as cr^ ^ = o"^ ^ — A. 



A previous study using data from the Swedish Vacuum Solar Telescope (Ivan Ballegooijen et al 



19981 ) assumed c(t) to be a Lorentzian. Here, we clearly see that c(t) differs from a Lorentzian, and 



^The correlation at zero-time lag is the variance of the velocity distribution. 
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it can be fitted with a function £, which is a sum of the true correlation of solar origin (C') and A, 
given by 

e:„(A,r,K) = e:;(T,K) + A„, (8) 



<(r,K) = a + — (9) 



where 



is a generalized Lorentzian. A, r (correlation time), and k (exponent) are the free parameters of 
the fit; a and b are the functions of (A,r, k), which are determined analytically by least square 
minimization (see Appendix|X]). We also bring to the notice of the reader that our formula for C is a 
monotonically decreasing function of t„. However, there is an unexplained increase in the observed 
Cyy beyond ±100 s (panel (b) in Figure [6]). To eliminate any spurious results due to this anomaly, 
we use a maximum delay time of ±105 s to fit c with <t by minimizing the sum of the squares of 
their difference, as defined in Equation (lAip . 



The top panel in Figure [7| shows the results listing the best fit values of the free parameters 
(A, T, k), a, and b for a maximum t„ of ±105 s. C (black) and C (thin red) are plotted as functions 
of the delay time over Cxx,n (left, symbols), and Cyy^n (right, symbols). The value of A where has 
its global minimum is found to be 0.75 km^ s~^, for both Cxx and Cyy. The bottom panel shows the 
contours of as a function of r and k at A = 0.75 km^ s~^, and the min{x^) is denoted by plus 
symbols. Dashed and solid lines are the regions of 1.5 and 2 times the minix^) respectively, is 
a well bounded function for k < 2, confirming a cusp-like profile. The correlation time is 22 — 30 
s, which is about 4 — 6 times the time cadence. 

Taking into account the variance in errors (i.e., A = 0.75 km^ s~^), we get dp = 3 km, and the 
corrected rms velocities (<t^,c) of Vx and Vy are now 1.00 and 0.86 km s~^. These results are plotted 
as dashed curves in panels (a) and (b) of Figure \5\ and the values are tabulated in the last row of 
Table [TJ The corrected distribution of v is shown as a dashed Rayleigh distribution in panel (c) . 
With higher cadence observations, these results can be refined and modified, as (A, r, k) depend on 
the shape of the core of c. By comparing the SST and Hinode results, we expect that the observed 
c probably increases rapidly below 5 s and thus changing the set of parameters to some extent. 



5. Summary and Discussion 



We studied the proper motions of the BPs using the wideband Ha observations from the SST 
and the G-band data from Hinode. BPs are manually selected and tracked using 5" x 5" areas 
surrounding them as reference frames. The quality of the SST observations allowed us to measure 
the BP positions to a sub-pixel accuracy with an uncert a inty o f only 3 km, which is at least seven 
times better than the value reported by iNisenson et al.l (120031). and co mparable to the rms value 



of 2.7 km due to image jittering reported by lAbramenko et al.l ()201ll ). They adopted this rms 



value of 2.7 km as a typical error of calculations of the BP position. We found that the horizontal 
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motions of the BPs in x and y are Gaussian distributions with raw (including the true signal and 
measurement errors) rms velocities of 1.32 and 1.22 km s~^, symmetric in u— space, observed at 5 s 
cadence. The above estimate of the measurement uncertainty is obtained from a detailed analysis 
of the velocity auto-correlation functions. For this, we fitted the observed c{t) with C, a function 
of the form shown in Equation ([8|), and estimated an rms error of about 0.87 km s~^ in and Vy. 
The removal of this error makes the Vx and Vy Gaussians narrower with new standard deviations 
1.00 and 0.86 km s~^ (a fractional change of 30%). The total rms velocity {vx and Vy combined) is 
1.32 km s^^. The correlation time is found to be in the range of 22 — 30 s. 

Following is a brief note and discussion on the additional results we derive from our work. 
BPs are advected by the photospheric flows. Thus, taking these features as tracers, we can derive 
the diffusion parameters of the plasma. As BPs usually have life times of the order of minutes, 
the motion of these features can be used to study the nature of photospheric diffusion at short 
time scales. The mean squared displacement of BPs ((Ar)^) as function of time is a measure of 
diffusion. It is suggested in the literature that ( (Ar)^) can be approximated as a power law with 



an index 7 (i.e., ((Ar)^) ^ , see for example ICadavid et al.lll999l : lAbramenko et al.ll201ll ). In 



Figure [HI we plot the observed ((Ar)^) (symbols), for the 200 s interval on a log-log scale. Solid 
lin e is the least squa r e fit w ith a slope of 1.59, which is consistent with the value 7 = 1.53 found 



bv lAbramenko et al.l (|201ll ) for quiet sun. Despite the differences in the observations (instruments 
and observed wavelengths), and analysis methods (identification and tracking of BPs), a close 
agreement in the independently estimated 7 suggests that this is a real solar signal. Both these 
results assert the presence of super-diffusion (i.e., 7 > 1), for time intervals less than 300 s. Since 
most of the BPs in this study are tracked for only 3 — 4 min, we cannot comment on the diffusion 
at longer times. 

Note that there is a general relationship between the mean squared displacement ((Ar)^) and 
the velocity auto-correlation function 

{(Arf) = {(^j vx{t')dt'^ ) + {(^j vy{t')dt'^ ) (10) 

t t 

= J ^'{t" -t')dt'dt", (11) 


where we assume isotropy of the BP motions {€-'xx = ^'yy = ^')- For a known auto-correlation or 
mean squared displacement, the other quantity can be derived using the above relation. 

We already saw that the horizontal motions of the BPs yield several important properties of 
the lower solar atmosphere. One more such important property is the possibility of the generation 
of Alfven waves due to these motions. Here we qualitatively estimate and compare the power 
spectrum of horizontal motions as a function of frequency for two forms of the velocity correlation 
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functior@: (a) the form C' (Equation dH)), obtained in this study, and (b) a Lorentzian function. 
For case (a) we use a = 0, and also assume that C^^.^^ = ^'yy^ni with the parameters b, r, and k 
taking the mean values of x and y. For case (b) we use a modified form of £' with k = 2. The other 
parameters (a, 6, and r) are the same as in case (a). Figure [9] shows the power spectra for the two 
described cases: (a) solid line, and (b) dashed line. We observe that for frequencies exceeding 0.02 
Hz (< 50 s), the horizontal motions generally have more power in case (a) as compared to case 
(b). This highlights the fact that the dynamics of the BPs at short time scales are very important. 
Therefore, it is highly desirable to do these observations and calculations at very high cadence. 

The measurements presented in this paper provide important constraints of models for Alfven 
and kink wave generation in solar magnetic flux tubes. As discussed in the Introduction, such waves 
may pl ay an important role in chroinospheric and coronal heating. In the A lfven wave turbulence 



model (Ivan Ballegooijen et al.l l201ll : lAsgari-Targhi &: van BallegooijenI l2012l ) it was assumed that 



the photospheric footpoints of the magnetic field lines are moved about with rms velocity of 1.5 
km s~^, similar to the rms velocity of 1.32 km s~^ found here. However, the models include only 
the internal motions of a flux tube, whereas the observations refer to the displacements of the flux 
tube as a whole. Clearly, to make more direct comparisons between models and observations will 
require imaging with high spatial resolution (< 0.1 arcsec). This may be possible in the future with 
the Advanced Technology Solar Telescope. 

In this work we presented the results of the BP motions, some of their implications and use 
in the context of photospheric diffusion and coronal wave heating mechanisms. We interpret the 
location of the intensity maximum of a BP as its position at any given time. This is certainly 
plausible for time periods when we begin to see the physical motion of a BP as a rigid body due to 
the action of the convection on the flux tubes. But at timescales shorter than one minute, other 
interpretations are also plausible: the motions marked by the intensity maxima could be intensity 
fluctuations in an otherwise static BP. Nevertheless, these fluctuations are manifestations of some 
disturbances inside the BP, which are equally important and interesting to explore further. 



^Fourier transform of the velocity auto-correlation is the power spectrum. 
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Table 1. Properties of the velocity distributions in Figure [5] 





Vx (km s ^) 


Vy (km s 




{Vx) 


(y{vx) 




a{vy) 


Histogram 


0.18 


1.58 


0.19 


1.54 


Gaussian Fit (ay^r) 


0.01 


1.32 


0.01 


1.22 


Corrected Distribution (c7i,,c) 


0.01 


1.00 


0.01 


0.86 
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A. Determination of a and h 



In this section, we briefly describe a method of determining a and h for a set of parameters 
(A,r, k). We define of autocorrelation functions of Vx and Vy as 



N 



xL(^>T,«^) = X] [cxx,n - ^n(A,r,K)]^and 

n=-Af 

X^j,(A,T,k) = [Cyy,n-^n{^,T,K)f (Al) 

n=-Ar 

where, Cxx,n and c^^^n are the observed autocorrelation values of velocities Vx and Vy and £ is a 
model of the correlation function given by Equation By minimizing the with respect to a 
and b (i.e., = and = 0, separately for x and y), and solving the resulting system of linear 
equations, we have 



where. 



a = 2n + 1 

N 



a' = /3 



„=-iv 1 + ( ¥ 



1 



n=-Af 



n=-N 
N 



n=-N 



1 + ( IM 
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2006-06-18 13:10:39 UT 




arcsec 



1. — First image from the time sequence of SST wideband Ha observations at 13:10 UT on 18 
2006. Black arrow is pointing towards solar north and the white arrow is towards disk center. 
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Fig. 2. — Illustration of different offsets seen in the full and partial FOV of Hinode data. Dashed 
curves in the left (right) panel show the drifts in the x— (y—) direction of the four selected quadrants. 
Thick dashed profile is the average of four dashed curves. Thick red profile is the drift of the full 
FOV (see text for details). 
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Fig. 3. — Examples of the paths of four BPs taken from the SST data. Initial position of each BP 
is marked with a star. Time shown at the top right corner in each panel is the duration for which 
respective BP is followed. 
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Fig. 4. — Velocities Vx and as a function of time for a typical BP (shown here for BP^3, see 
Figure [3] for the track of BP^^S). 
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Fig. 5. — Histograms of measured BP velocities: (a) Vx, (b) Vx, (c) f = y + v"^. Solid Gaussians in 
the top panels are fits to the histograms. Dashed Gaussians in panels (a), (b), and dashed Rayleigh 
profile in panel (c), are the new distributions of velocities after correcting for the measurement 
errors (see text for details); (d) shows Vx plotted against Vy. 
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Fig. 6. — Correlation functions of BP velocity Vx and Vy. (a) Observed auto-correlation Cxx,n as a 
function of delay time t (black: SST, red: Hinode). (b) Similar for the observed auto-correlation 
Cyy,n- (c) Cross-correlatiou Cxy,n as function of delay time, (d) Number of measurements per bin 
used in panels (a), (b), and (c). 
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Fig. 7. — Top: £ (black curve) plotted as a function of delay time with the best fit values of a, b, 
A, r, and k obtained by minimizing the (see Appendix|A]), of the observed c (shown as symbols, 
Cxx- left, and Cyy: right, also shown as black curves in the top panels of Figure [6|), and the modeled 
correlation function (t, for a delay time of ±105 s in steps of 5 s. Thin red curve is the profile of 
Bottom: Contour plots of a function of k and r, for a value of A (A^; = A^ = 0.75 km^ 

s~^), where attains the global minimum. Plus symbol is the global minimum of x^; dashed and 
solid lines are the contours of 1.5min{x^) and 2min{x^) respectively. 
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Fig. 8. — Mean squared displacement ((Ar)^) as a function of time t on a log-log scale. Solid line 
is the least square fit of the observations (symbols), with a slope of 1.59. 
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Fig. 9. — Power spectrum of the horizontal motions of BPs as a function of frequency derived from 
auto-correlation function for two cases. Solid line: case (a) - from this study. Dashed line: case (b) 
- from a Lorentz profile with same a, b, and r as in case (a) but with k = 2 (see text for details). 



